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Abstract The current flow along the boundary of graphene stripes in a perpendicular magnetic field is 
studied theoretically by the nonequilibrium Green’s function method. In the case of specular reflections 
at the boundary, the Hall resistance shows equidistant peaks, which are due to classical cyclotron motion. 

When the strength of the magnetic field is increased, anomalous resistance oscillations are observed, similar 
to those found in a nonrelativistic 2D electron gas [New. J. Phys. 15:113047 (2013)]. Using a simplified 
model, which allows to solve the Dirac equation analytically, the oscillations are explained by the inter¬ 
ference between the occupied edge states causing beatings in the Hall resistance. A rule of thumb is given 
for the experimental observability. Eurthermore, the local current flow in graphene is affected significantly 
by the boundary geometry. A finite edge current flows on armchair edges, while the current on zigzag 
edges vanishes completely. The quantum Hall staircase can be observed in the case of diffusive boundary 
scattering. The number of spatially separated edge channels in the local current equals the number of 
occupied Landau levels. The edge channels in the local density of states are smeared out but can be made 
visible if only a subset of the carbon atoms is taken into account. 


1 Introduction 

Nowadays, graphene is maybe the most studied material 
in condensed matter physics because of its numerous ex¬ 
ceptional properties and their potential technological ap¬ 
plications, see w and references therein for an overview. 
In particular the electronic transport of charge carriers in 
graphene is of enormous interest due to the promise of 
novel electronic devices like foldable displays and high- 
frequency transistors [5]. It has also been shown recently 
that strain and deformation of a graphene stripe, as caused 
by the absorption of atoms for example [7], give rise to a 
strong pseudo-magnetic field HIE], which affects signifi¬ 
cantly the current flow m- 

In this paper, we study magnetotransport along the 
edges of graphene stripes. As sketched in figure elec¬ 
trons are injected coherently at one point S on the bound¬ 
ary of the graphene stripe and focussed by a homogeneous 
perpendicular magnetic field B onto another point Pi on 
that boundary. In the classical regime (blue trajectories) 
resonances are observed, if a multiple of the cyclotron di¬ 
ameter equals the distance between the injecting and col¬ 
lecting point contacts mma. For large Fermi wavelength 
and long phase coherence length, additional interference 
effects appear. This regime of coherent electron focusing 
has been studied for the first time by van Houten et al. 
in a nonrelativistic two-dimensional electron gas (2DEG) 
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m but it has become a topic of current interest again, 
since the first focusing experiments in graphene have be¬ 
come possible [Hill]. The magnetic focusing in graphene 
pn junctions has been studied theoretically m and snake 
states at such a pn interface have been predicted [IIlIIH]. 
It has also been suggested to study by coherent electron 
focusing the structure of graphene edges HSj. Recently, 
also the effects of disorder [20] and spin-orbit interaction 
[2IH25] have been investigated in a nonrelativistic 2DEG. 
On the other hand, graphene stripes in a strong magnetic 
field show the quantum Hall effect EgEii, which is ex¬ 
plained by the transport through edge channels along the 
boundary of the system, see the red lines in figure 

Here, we study theoretically the system properties from 
the classical to the quantum regime. In particular, we dis¬ 
cuss the novel effects which emerge, when the two regimes 
are bridged by suitable system parameters. If the scatter¬ 
ing at the boundaries is specular, we observe in this inter¬ 
mediate regime anomalous resistance oscillations, which 
are neither periodic in B (classical cyclotron motion) nor 
periodic in 1/5 (quantum Hall effect). Using a simplified 
model, which allows to solve the Dirac equation analyti¬ 
cally, we explain these oscillations by the interference of 
the occupied edge channels. These anomalous resistance 
oscillations have been reported recently in a nonrelativis¬ 
tic 2DEG [28]. Beyond this, in graphene the local current 
flow is affected significantly by the boundary geometry. 
We show that on armchair edges a finite edge current is 
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assume that the influence of the temperature is negligible 
and thus, set it to zero. We consider graphene stripes with 
a zigzag boundary in between S and Pi, as well as with an 
armchair boundary, see the inset of figure Other pos¬ 
sible edge reconstructions, see e.g. |30]-[32], are not con¬ 
sidered here. Note that the orientation of the graphene 
lattice with respect to the used coordinate system is not 
changed in the two stripes. Armchair edges run along the 
x-axis, whereas zigzag edges are oriented along the ^-axis, 
see the inset of figure]^ 


Figure 1. Magnetotransport in graphene stripes is stud¬ 
ied. Electrons are injected at the source S and focussed by 
a perpendicular magnetic field B onto the contact Pi. We 
calculate the current Isd between source S and drain D as 
well as the voltage drop Up-^p^ between the voltage probes 
Pi and P 2 . Using these quantities, we study the gener¬ 
alized Hall resistance R^y = Bp^p^jlsD as a function of 
B. Cyclotron orbits at low magnetic field are sketched by 
the blue trajectories. The edge channel transport of the 
quantum Hall effect at high magnetic field is indicated by 
red lines. 


3 Calculations 

In this section, we begin with a short introduction into the 
nonequilibrium Green’s function (NEGF) method, which 
is applied to study quantitatively the magnetotransport in 
graphene stripes. Afterwards, we solve the Dirac equation 
in a magnetic field for a graphene sheet bounded by a sin¬ 
gle infinite potential wall. This simplified model will help 
us to get insight into the results of the Green’s function 
calculations. 


present, while on zigzag edges the current is shifted to the 
interior of the stripe and vanishes exactly on the edge. 
We also give a rule of thumb for the experimental ob¬ 
servability of these effects. The quantum Hall staircase is 
observed in the case of diffusive scattering at the bound¬ 
aries, which guarantees that the phase coherence length is 
shorter than the relevant geometric lengths (point contact 
distance, system size). We show that the number of spa¬ 
tially separated edge channels in the local current equals 
the number of occupied Landau levels. In the local density 
of states (LDOS) the edge channels are smeared out but 
can be made visible if only a subset of the carbon atoms 
is considered. 


2 System 


We study a graphene stripe with a size of 140 nm x 90 nm, 
see figure [2 Metallic contacts with a width of 3 nm are at¬ 
tached at the edges of the stripe separated by a distance 
of P = 110 nm (measured between the middle of the con- 
tacs). We calculate the current Isd flowing between source 
S and drain D due to an infinitesimal bias voltage, as 
well as the voltage drop Up^p^ between the voltage probes 
Pi and P 2 . Using these quantities, we study the gener¬ 
alized Hall resistance Rxy = Up^p^lI sd as a function of 
an homogeneous perpendicular magnetic field B = —Bcz- 
The Fermi energy is set to /i = 0.06t = 168 meV, where 
t = 2.8 eV is the coupling between nearest neighboring 
carbon atoms at a distance a = 0.142 nm [3]. In the exper¬ 
iment, usually the electron density is constant while the 
chemical potential is oscillating. However, this would only 
slightly displace the transitions between the Hall plateaus 
and would not qualitatively change our results, see also 
m- For simplicity, we do not take into account the Zee- 
man spin splitting, see remarks in section 4.4 We also 


3.1 The nonequilibrium Green’s function method 

We start from the tight-binding Hamiltonian of the gra¬ 
phene stripe 


( 1 ) 

(d) 

where (ij) means nearest neighbors at a distance a with 

coupling t. The are the Pz orbitals of the carbon 

atoms on sublattice A and B, respectively, see the red and 
black marked atoms in the inset of figure In graphene 
nanoribbons it can be necessary to take into account also 
the interaction to second and third nearest neighbors 133 - 
[ 35 ] . However, for the graphene stripes studied here, it is 
sufficient to consider only nearest neighbors, as our main 
results remain qualitatively unchanged if also second and 
third nearest neighbors are taken into account. 

In the tight-binding Hamiltonian the effect of the mag¬ 
netic field B is taken into account by the Peierls substi¬ 
tution [36] 

= ( 2 ) 

where A is the vector potential of the magnetic field. The 
path integral is along the straight connection between the 
position of carbon atom i and j. 

The Green’s function of the graphene stripe is defined 
as |37H39] 


G = 


E-H-EtiD 


( 3 ) 


where E is the energy of the charge carriers. The influence 
of each of the Nc contacts is taken into account by an 
imaginary self-energy 

Vi 


( 4 ) 
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with broadening 77 = 1.25t = 3.5 eV, representing metallic 
contact regions. The sum is over all carbon atoms which 
are coupled to the same contact k. 

The transmission from contact j to contact i is then 
given by 


Tij = 4Tr (Im (/^) GIm (G,-) G^) . (5) 

and the total current at the ith contact reads 

~ ~ mO 5 (^) 

j 

where is the chemical potential of contact i and j, 
respectively. The generalized Hall resistance is then given 
by 


_ upiP2 _ (^Pij ~ '^P2j) Tjs 

IsD 2e^ Tds + Yhij TDiTlijTjs ’ 


where 


The sums in 0 are over the contacts with unknown chem¬ 
ical potential, whereas the sum in ^ is over all contacts 
including source and drain. 

The local current of electrons, which originate from the 
source with energy /i and which flow from atom j to the 
neighboring atom i, is given by iiniiii] 


—T- 


E 


k^i 


Tik 


i = j. 


(8) 


2p 

hj = -Im {h;,A%) , (9) 

where the Hij are the matrix elements of the Hamilto¬ 
nian Q . The spectral function for electrons from the source 
is defined as 


=--GIm(rs)G+. ( 10 ) 

TT 

The diagonal elements of the spectral function give the 
local density of states (LDOS), which is accessible to these 
electrons. 

Finite system size effects, such as standing waves be¬ 
tween the system boundaries, would distort the magne¬ 
totransport strongly. Therefore, diffusive boundaries are 
used at those edges, which are not important for the focus¬ 
ing experiment, see the dashed edges in figureDiffusive 
boundaries are implemented mathematically by additional 
virtual voltage probes, which randomize phase and mo¬ 
mentum of the charge carriers and thus, suppress standing 
waves in the system. The chemical potential of the virtual 
reservoirs is determined by the condition that no charge 
carriers can be gained or lost at a virtual reservoir (current 
conservation constraint), see [28l|42j[43] for details. 

3.2 Dirac equation in a magnetic field 

For Fermi energies close to = 0, the physics of graphene 
takes place at two points K and K' in momentum space. 


At these points the dispersion relation is linear E{k) = 
hvp \k\ and the charge carriers behave as relativistic mass¬ 
less particles described by the Dirac Hamiltonian 


Hkik' = vf 


0 Px T '^Py 
Px ± ipy 0 


( 11 ) 


with vf = 3atl2h. Note that in order to derive from the 
tight-binding Hamiltonia n ^ the Dirac Hamiltonian (11) 
in its common notation [TJ [SflS] , we also applied an unitary 
transformation 


U = 


1 0 
0 i 


( 12 ) 


The effect of the magnetic field is taken into account 
by minimal gauge invariant coupling p ^ p — eA, where 
A = BySx is the vector potential chosen for the armchair 
stripe, see figure]^ (bottom). To solve the Dirac equation, 
we insert the two linear equations into each other, keeping 
in mind that [px^Py] = —leBh. By means of the ansatz 
'0 b(^) = we obtain the Schrodinger equation 

of a harmonic oscillator 


ExB{y) = 


—+ {y-Vk) 


XB{y), (13) 


which is shifted by yk = and rescaled in energy E = 
2 ^ + ^ with ojc= ^ and = ^ - Thus, the eigenen- 
ergies are given by 


Ej, = ± 


yj2mvp hujc 


z/ > 0. 


(14) 


The eigenstates on sublattice B read xb ^ ), 

where Vy{ii) = Dy{x^y )/\/zd are rescaled parabolic cylin¬ 
der functions [44]. The eigenstates on sublattice A follow 
directly from xb^ the Dirac equation and the recursion 
relation {dy + y) Vj^{y) = x^Vjy-i{y), see [44|. The solu¬ 
tion of the Dirac equation at the K' valley can be obtained 
easily by interchanging the two sublattices, see (11). The 
eigenenergy spectrum is unchanged and hence, twofold de¬ 
generate. Thus, the eigenfunctions are given by 


where is a normalization constant, ^ = [y — /£b 

and Vx<o = 0. The different signs of the eigenstates cor¬ 
respond to the signs of the eigenenergies. We also applied 
the unitary transformation ( [T^ in order to get the correct 
phase between the wavefunctions on the sublattices. For 
the zigzag stripe, see figure]^ (top), we choose the vector 
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potential A = —BxCy^ to get the eigenstates 

where ( = (x + /c) • In an infinitely extended sys¬ 

tem, the index u has to be an integer n = 0,1, 2 ... (Lan¬ 
dau level index), because the eigenfunctions have to be 
normalizable. In this case, the parabolic cylinder functions 
can be simplified by Vn{y) = e“^ j where 

Hn{y) are the Hermite polynomials. 


3.3 Dirac equation with an edge in a magnetic field 

To understand the magnetotransport in graphene stripes, 
we solve the Dirac equation bounded by an edge under 
the effect of a magnetic field. In general, the solution of 
the Dirac equation is given by a linear combination of the 
solutions at both valleys 


= Cl e'^"^tpK{r) + C2 (17) 

where ci and C 2 are complex constants. At the zigzag edge, 
we obtain by means of 


^ 7 


,(r) = 




+C2e 


±r>^-i(C) 


(18) 


As only carbon atoms of one sublattice appear at a zigzag 
edge, see the inset of figure (top), the wave function has 
to vanish only on one of the two sublattices. The condition 
= 0) =0 leads to the two solutions 

ci = l, C 2 = 0: = 0, (19a) 

Cl =0, C2 = 1 : V^{£Bk) = 0. (19b) 


Thus, for given iBk = ^hk‘^/eB the index y is deter¬ 
mined by the zeros of the rescaled parabolic cylinder func¬ 
tions. The first set of solutions is located at the K valley, 
whereas the second set is located at the K' valley. The 
resulting energy bands (14) are depicted in figure(top). 
At large k the discrete Landau levels for integer values 
of y = n can be observed. The distance of these Landau 
levels decreases with ^/n. When the apex of the parabola 
Vk = £%k approaches the wall by decreasing k, the energy 
bands are bent upwards and their degeneracy is lifted. 
Also a dispersionless state Ej^=o = 0 at the K valley can 
be seen. The occupied edge states at the Fermi energy 
(dashed horizontal line) are indicated by dots. 


At an armchair edge, we obtain by means of (15) 


S'ac(r) = 




Cl e 


+C 2 e 




T^AO 


( 20 ) 




Figure 2. Energy bands of graphene bounded by an edge 
in a magnetic field of B = 15.6 T. At large k we observe 
discrete Landau levels. When k is decreased, the energy 
bands are bent upwards and their degeneracy is lifted. In 
the case of a zigzag edge (top) the red curve indicates 
solutions at the K valley, while the blue curve gives so¬ 
lutions at the K' valley. In the case of an armchair edge 
(bottom) the valleys are mixed, which leads to two sets of 
solutions indicated by the blue and red curve. The occu¬ 
pied edge states at the Fermi energy (dashed horizontal 
line) are marked by blue and red dots. 


As at armchair edges both sublattices appear, see the inset 
of figure]^ (bottom), the wave function has to vanish on 
both of them. The condition EA{y = 0) = EB{y = 0) = 
0 requieres that the coefficient determinant of the linear 
equation system for ci and C 2 vanishes 

V,_i{-iBk) T V,{-iBk) = 0 (21) 

and leads to the solutions 

Cl = 1, C2 = ±1. (22) 

Thus, at an armchair edge both valleys are intermixed. 
The two eigenenergy bands in figure (bottom) show not 
only that their degeneracy is lifted in vicinity of the edge 
but also shallow valleys, which are not present at a zigzag 
edge. The solution of the Dirac equation at zigzag and 
armchair edges in a magnetic field can also be found in 

HSHIH]. 
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4 Results and discussion 

4.1 Density of states 


Let us discuss briefly the density of states (DOS) in the 
studied devices, which is shown in figure For energies 
E > 0.02t the DOS in the zigzag and the armchair stripe 
agree well with the DOS Ddi(L^) = from the Dirac 

Hamiltonian. However, the zigzag stripe shows a distinct 
peak at = 0, which cannot be observed in the case of 
an armchair stripe. This peak can be attributed to the 
dispersionless state shown in figure]^ (top), which does 
not contribute to electron transport and is located on the 
surface of the stripe. A surface state is possible at a zigzag 
edge, because only carbon atoms of a single sublattice ap¬ 
pear there. Thus, at the edge the wave function has to 
vanish only on one sublattice, while the surface state re¬ 
sides on the other sublattice. At an armchair edge atoms 
from both sublattices appear and a surface state is not 
possible, see 0 EH ESI EOj for details. However, in fig¬ 
ure the DOS of the armchair stripe is also nonzero at 
£■ = 0. These are states induced by the contacts EH, 
which contribute to the observed finite conductivity of 
graphene at the Dirac points [26l |27l [52H54] . As a con¬ 
sequence of this, at /i = 0.06t = 168 meV, the carrier 
densities in the zigzag stripe Uzz = 3.7 • 10^^ cm“^ and 
the armchair stripe rzac = 3.3 • 10^^ cm“^ are somewhat 
higher than expected from the linear DOS of the Dirac 
Hamiltonian n^i = 2.5 • 10^^ cm“^. 



Figure 3. Average DOS in the studied graphene stripes. 
The blue curve gives the DOS in the zigzag stripe, while 
the red curve gives the DOS in the armchair stripe, see the 
inset of figure The DOS of both stripes agrees well with 
the DOS from the Dirac Hamiltonian (green curve). The 
peak in the DOS of the zigzag stripe at = 0 is caused 
by the dispersionless surface state, which does not exist in 
armchair stripes. However, also the armchair stripe has a 
nonzero DOS at = 0, because of contact induced states. 
The dashed vertical line indicates the Fermi energy. 
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Figure 4. Hall resistance Rxy as a function of the mag¬ 
netic field B for the zigzag stripe (top) and the armchair 
stripe (bottom). The blue curve gives Rxy in the case of 
specular reflections at the boundary between S and Pi, 
whereas for the red curve the scattering at this bound¬ 
ary is diffusive. The Hall resistance starts with peaks, 
which can be understood by classical cyclotron orbits (23), 
see the dashed vertical lines. In a strong magnetic field 
B > 10 T, we observe superimposed upon the quantum 
Hall plateaus anomalous resistance oscillation, which can¬ 
not be explained by cyclotron orbits. 


4.2 Cyclotron motion and the quantum Hall effect 

The Hall resistance Rxy as a function of the magnetic field 
P, calculated by means of the NEGF method ([^, is shown 
in figure In the case of specular reflections at the bound¬ 
ary between S and Pi (blue curve), the Hall resistance of 
both stripes shows at low magnetic field 0T<P<10T 
a series of equidistant peaks located approximately at 

Bn = ^^n, n = 1,2,3,..., (23) 

evpL 

see the dashed vertical lines. At these magnetic fields a 
multiple n of the cyclotron diameter 2 |p| /eB equals the 
distance L between injector and collector. Cyclotron or¬ 
bits can be clearly seen in figure which shows the local 
current and the local density of states (LDOS) of electrons 
originating from S with energy ji. Note that the shown lo¬ 
cal current and the LDOS have been averaged over the 
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Figure 5. Local current (arrows) and LDOS (shading) of electrons originating from S with energy fi. In the zigzag 
stripe (left column) and the armchair stripe (right column) cyclotron orbits can be clearly seen. At the armchair edge 
a distinct edge current can be observed, which is not present at the zigzag edge. Note that the shown local current 
and the LDOS have been averaged over the honeycomb cells. 
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Figure 6. The edge channel transport of the quantum Hall effect can be observed clearly, when the direction of the 
magnetic field is reversed. A finite current flows at the armchair edge (right), whereas the current vanishes at the 
zigzag edge (left). This can also be seen in the transverse current through the dashed vertical line in figureIn the 
LDOS only a single broadened edge channel can be recognized, instead of spatially separated edge channels. 


honeycomb cells. These current flow paths can be mea¬ 
sured by scanning tunneling microscopy [55] . 


Extended quantum Hall plateaus in a strong magnetic 
field H > lOT can be observed, if the boundary in be¬ 
tween S and Pi is diffusive (red curve), or if the magnetic 
field is reversed and the current passes by th e other diffu¬ 
sive boundaries (see the end of section 3.1). The current 
is carried through edge channels along the boundaries, see 
figure The quantum Hall effect in graphene can be un¬ 
derstood easily by the eigenenergy spectra shown in fig¬ 
ure The number of the occupied edge states at the Fermi 
energy equals 2n+l, where n is the Landau level index. As 


every occupied edge state is a ballistic conductor, which 
contributes with 2e^ jh to the total conductance, the Hall 
resistance reads Rxy = ^ 2n-\-i ^ ^ ~ 0,1,2,.... This ex¬ 
plains the quantum Hall staircase observed in figure 
which is one of the definitive fingerprints of a relativis¬ 
tic 2DEG [26l [23 [56ll58] , because it differs significantly 
from the nonrelativistic case In figure we can also 
observe that the transitions between the Hall plateaus dif¬ 
fer slightly in the two stripes. This can be explained by 
the shallow valleys in the band structure at an armchair 
edge, which are not present at a zigzag edge or when the 
scattering at all boundaries is diffusive. 
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4.3 Anomalous resistance oscillations 

In the case of specular scattering between S and Pi and 
in magnetic fields B > 10 T, we observe - superimposed 
upon the quantum Hall plateaus - anomalous resistance 
oscillations, which cannot be understood by classical cy¬ 
clotron motion. In particular, when only two Landau level 
are occupied (16 T < B < 26 T), the oscillations become 
very clear and regular. Their frequency increases rapidly 
whenever a Landau level is pushed towards the Fermi en¬ 
ergy and a transition between Hall plateaus appears (com¬ 
pare blue and red curves in figure]^. Finally, the oscilla¬ 
tions vanish completely, when only a single edge channel 
is occupied {B > 26 T), and the Hall plateau Rxy = 1 ap¬ 
pears (not shown in figure [^. These resistance oscillations 
can be understood by means of the solution of the Dirac 
equation. In the zigzag stripe the edge states are given by 
(18) and (H^. We superimpose the plane wave part of the 


occupied edge states 
IV-grl^ = - 


E n+1 

i=l 




S', Pi 



{2n + iy 

(24) 

where {')s,Pi means spatial averaging over the finite width 
of the injector and collector contacts. The occupied edge 
states in the K valley are denoted by ki and the states 
in the K' valley by qi, see the red and blue dots in fig¬ 
ure The normalized absolute square of these superim¬ 
posed plane waves agrees almost perfectly with the NEGF 
calculation, see figure]^ (top). Thus, all focusing peaks can 
be understood by the interference of the plane wave part 
of the occupied edge states. The anomalous resistance os¬ 
cillations are beatings, which appear when only some few 
edge channels are occupied. These beatings are very clear 
and regular if only two Landau levels are occupied. Their 
frequency increases rapidly, whenever the highest occu¬ 
pied Landau level approaches the Fermi energy, because 
its intersection point with the Fermi energy and thus, the 
corresponding /Cmax (oi* ^max) increases strongly. The dif¬ 
ference of /Cmax to the other, much smaller kn leads to a 
high frequency beating. Finally, when only a single edge 
channel is occupied, the beating and thus, the oscillations 
in the Hall resistance vanish. In the armchair stripe, the 
solution of the Dirac equation is more complicated, see 
(20), ( [^ , and ( [^ , because the valleys are intermixed. 
We found best agreement to our Green’s function calcula¬ 
tions, see figure [^(bottom), if we use also for the armchair 
stripe where the ki and Qi denote the two sets of so¬ 
lutions. 

In a magnetic field B < 16 T, the simplified model 
shows smaller highly oscillating peaks, which are due to 
the interference of numerous plane waves. These highly 
oscillating peaks are more pronounced in the armchair 
stripe, where additional interference between the K and 
K' valley takes place, which is not present in the zigzag 
stripe, compare and ( p^ . The highly oscillating peaks 
are not present in the NEGT calculations due to the diffu¬ 
sive boundaries. Probably they neither appear in the ex¬ 
periment due to the presence of decoherence. In a stronger 



Figure 7. Normalized absolute square of the superim¬ 
posed plane wave part of the occupied edge channels (24), 
(red curve) agrees well with the Green’s function calcu¬ 
lation of the Hall resistance (blue curve). All resistance 
oscillations can be understood by the interference of the 
edge channels. The anomalous oscillations are beatings, 
which appear when only some few edge channels are oc¬ 
cupied. 


magnetic field, the highly oscillating peaks disappear and 
the simplified model agrees very well with the NEGF cal¬ 
culations, because the superposition of only few eigen¬ 
states (see figure]^ leads to beatings. The oscillation are 
almost independent from the edge geometry, apart from 
slight differences in their frequency and phase. 

The beatings, which appear in the case of only two oc¬ 
cupied Landau levels, can be used to determine precisely 
the distance between the injector S and collector Pi. In 
figure the almost perfect match of the positions of all 
extrema in the range 13T < H < 26T is obtained only, 
if L = 110 nm is chosen in (24) for the distance between 
S and Pi. In order to explain, why in armchair stripes 


the classical focusing peaks deviate slightly from their ex¬ 
pected positions, see figure (bottom), we could assume 
hypothetically a slightly larger distance L = 120 nm be¬ 
tween injector and collector. In this case, the classical fo¬ 
cusing peaks would appear exactly at the expected posi¬ 
tions, but the beatings would absolutely not fit to (24). 
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Figure 8. Hall resistance of larger graphene stripes calculated by means of (24). Contacts with a width of 20 nm are 
attached at a distance of 450 nm. Right: When the Fermi energy is set to /i = 260 meV, corresponding to a carrier 
density of rigr = 6.0 


10^^ cm 


only classical equidistant focusing peaks can be observed. As reported in [19], the 
classical focusing peaks of higher order (n > 4) are clearly visible at armchair edges (top) but are suppressed at zigzag 
edges (bottom). Left: When Fermi energy /i = 80meV and carrier density rigr = 5.7 • 10^^ cm“^ are lowered, we find 
classical equidistant focusing peaks followed by anomalous oscillations. 


Also finite size effects can be ruled out as these deviations 
are not present in zigzag stripes of the same size. One rea¬ 
son for the shift of the classical focusing peaks could be 
the distinct edge current observed only at armchair edges 
or edge dependent scattering m- 

Although the charge carriers in graphene behave as 
relativistic massless fermions, the studied stripes show 
properties similar to a nonrelativistic 2DEG [28]: Clas¬ 
sical focusing peaks in weak magnetic fields, followed by 
anomalous resistance oscillations when the magnetic field 
strength is increased. In both systems the resistance os¬ 
cillations can be explained by the interference of the plane 
wave part of the occupied edge states. However, in graphene 
the linear dispersion, the valley degeneracy (symmetry 
points K and K' in momentum space) as well as the non¬ 
trivial edge geometry add subtle but important new as¬ 
pects. In this way, at first sight the local current flow looks 
similar in both systems (cyclotron orbits, edge channels), 
compare figure|^with figure 2 in [28]. However, the bound¬ 
ary geometry has a distinct effect on the local current flow, 
see figure which will be discussed in section [^ 


4.4 Experimental observability 

Due to computational limitations, the studied stripes are 
relatively small (L = 110 nm) and the considered mag¬ 
netic fields are quite strong (Hmax = 30 T). In these strong 
fields, also the Zeeman spin splitting of the Landau lev¬ 
els can be relevant [58] (60] [61] but we do not expect that 
the spin splitting changes qualitatively our findings. Al¬ 
though it is technically possible to realize such system 
parameters, this is not essential to observe our findings in 
an experiment. The important factor in an experiment is 
the maximal number of resolvable focusing peaks nmax, 
which is limited due to decoherence and partial diffusive 


scattering at the boundary. In order to observe anoma¬ 
lous resistance oscillations due to the interference of some 
few edge channels, the distance L between injector and 
collector as well the Fermi energy /i have to be tuned in 
such a way that the maximal number of possible specular 
reflections fulfills the rule of thumb 


1 L II 
6 a 


(25) 


which can be derived easily by (14) and (23). Of course, 
mean free path and phase coherence length also have to 
be comparable with L. To our knowledge in most focusing 
experiments such system parameters have been used that 
the regime of coherent electron focusing and the quan¬ 
tum Hall effect are well separated, see e.g. Figure 10 in 
[13]. However, signs of the anomalous oscillations can be 
observed in different geometries [62] [63] . 

Because of the excellent agreement of the simplified 
model ( |2^ and the Green’s function calculations, see fig¬ 
ure we can use this simplified model to study larger 
stripes, for which NEGF calculations are demanding. We 
consider stripes at which 20 nm wide contacts are attached 
at a distance of 450 nm. This is approximately the same ge¬ 
ometry used in the recent focusing experiment in graphene 
M as well as in a theoretical study m- When the Fermi 
energy is set to /i = 260 meV corresponding to a car¬ 
rier density of ngr = 6.0 • 10^^ cm“^, the system is in the 
regime of classical equidistant focusing peaks (rimax ^ 49), 
see figure (right). In agreement with results reported 
by Rakyta et al. [19], the focusing peaks of higher order 
(n > 4) are clearly visible at armchair edges but are sup¬ 
pressed at zigzag edges. When Fermi energy fi = SO meV 
and carrier density Ugr = 5.7 • 10^^ cm“^ are lowered, we 
bridge the regime of coherent electron focusing and the 
quantum Hall regime (umax ^ 15), see figure M (left). 
The Hall resistance starts with equidistant classical peaks, 
but anomalous oscillations follow when the strength of the 
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Figure 9. Transverse current (blue curves) through the 
dashed lines in figure A finite current flows on the arm¬ 
chair edge (bottom) but the current vanishes on the zigzag 
edge (top). The current, calculated by the eigenstates of 
the Dirac equation (red curve), agrees with the NEGF 
calculation and allows to attribute the different edge cur¬ 
rents to the different boundary conditions of the stripes. 
We can identify two spatially separated edge channels, 
which equals the number of occupied Landau levels (with 
E > 0). Due to the boundary conditions, the edge chan¬ 
nels are more densely packed in the zigzag stripe. 

magnetic field is increased. This gives us confidence that 
the predicted resistance oscillations can be observed ex¬ 
perimentally. 


Figure 10. Energy resolved transverse current through 
the dashed lines in figure Warm colors indicate a cur¬ 
rent from Pi to S while cold colors correspond to a current 
in the opposite direction. As shown by the shading close 
to the edges of the stripes (i.e. close to x = 0 and ^ = 0, 
respectively), a finite current flows on the armchair edge 
(bottom), which is not present on the zigzag edge (top). 
The number of spatially separated edge channels equals 
the number of occupied Landau levels (with E > 0), al¬ 
though the two edge channels closest to a zigzag edge are 
hardly distinguishable. Surprisingly, close to the Landau 
levels regions of counterpropagating current can be ob¬ 
served (blue regions). However, the total (integrated) cur¬ 
rent is quantized and does not change its sign. 


4.5 Edge current flow 

In figures and we observe that a finite current flows 
on the armchair edge, whereas the current vanishes on 
the zigzag edge. This can be seen clearly in figure (blue 
curve), which shows the transverse current through the 
dashed vertical lines in figure It can be understood, if 
we calculate the transverse current by means of the eigen¬ 
states of the Dirac equation 

2n+l 2n+l 

/di(r) (X E 

i=l i=l 


where q is a normalization constant and the sum is over 
the occupied edge states, see the dots in figureAt zigzag 
edges the parabolic cylinder functions have to be zero, see 
(19), which results in zero edge current. At armchair edges, 
the sum of the parabolic cylinder functions has to be zero, 
see (21), which allows for a finite edge current. The trans¬ 
verse current calculated by means of (26) agrees well with 
the Green’s functions calculations, see the red curves in 
figure In the transverse current at H = —15.6T, we 
can identify two spatially separated edge channels. Thus, 
the number of spatially separated edge channels in the lo¬ 
cal current, averaged over the honeycomb cells, equals the 
number of occupied Landau levels (with energy E > 0), 
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compare with figure The lifting of their degeneracy at 
the edge is not resolved in the local current. Due to the 
boundary conditions, in the zigzag stripe the two edge 
channels are more densely packed and harder to separate 
than in the armchair stripe. Note that the total edge cur¬ 
rent is approximately independent from the edge geom¬ 
etry. The energy resolved transverse current in figure pT] 
confirms these findings. Surprisingly, it also shows coun- 
terpropagating currents close to the Landau levels, see the 
blue shaded regions, in which the current flows in the op¬ 
posite direction as in the red shaded regions. However, 
note that the total (integrated) current is quantized and 
does not change its sign. The counterpropagating currents 
are also found when the NEGF method is applied to a 
nonrelativistic 2DEG. At this point, their origin is not un¬ 
derstood, but they are also observed by Wang et al. [48] 
using the eigenstates of the Dirac equation. Also the de¬ 
pendency of the current on the edge geometry is reported 
in their work. Beyond that, we show in figure that a 
distinct armchair edge current appears also in the regime 
of coherent electron focusing. This distinct armchair edge 
current in focusing experiments could be measured exper¬ 
imentally by means of an additional voltage probe placed 
on the stripe’s edge or by contacting edge channels indi¬ 
vidually, as in [661167]. 


4.6 Local density of states 

The local density of states (LDOS) in figure averaged 
over the six carbon atoms of the honeycomb cells, shows 
only a single broadened edge channel. This can be seen 
clearly in figure]^ (black curves), which gives the LDOS 
along the dashed vertical line in figure In order to make 
individual edge channels visible in the LDOS, we have to 
select only a subset of the carbon atoms, see the blue and 
red curves for which only the atoms marked in the inset 
are taken into account. Note that in the armchair stripe 
two subsets give numerically identical results, see the blue 
curve. When every carbon atom is considered individu¬ 
ally, the LDOS oscillates rapidly between the blue and 
red curves in figure These oscillations have been re¬ 
ported in theoretical studies [45] |64| |68l [69], but to our 
knowledge an experimental confirmation is missing. The 
energy resolved LDOS, calculated numerically by means of 
the NEGF method, is depicted in figure [T^ Far from the 
edge the discrete Landau levels can be observed clearly. 
If the LDOS is averaged over the honeycomb cells (left 
column), the bending of the energy bands can hardly be 
discerned. It becomes more visible, if only a subset of the 
carbon atoms is taken into account (middle and right col¬ 
umn). In this way, we can observe how in the zigzag stripe 
(top row) the zeroth Landau level at = 0 splits into a 
dispersive edge state on the sublattice B (right) and a non- 
dispersive surface state on the sublattice A (middle). This 
surface state is not present in the armchair stripe. Similar 
results can also be obtained by means of the eigenstates 
of the Dirac equation, see m- Anyway, in the experiment 
it is not possible to select a subset of the carbon atoms. 


Thus, the measured LDOS looks similar to the figures in 
the left column, see |70]. 




Figure 11. LDOS of the zigzag stripe (top) and the arm¬ 
chair stripe (bottom) along the dashed vertical lines in 
figure The LDOS averaged over the six carbon atoms 
of the honeycomb cells (black curve) shows a single broad¬ 
ened edge channel. Individual edge channels become visi¬ 
ble, when only a subset of the atoms is taken into account, 
see marked atoms in the legend. Note that in the armchair 
stripe two subsets give numerically identical results (blue 
curve). 


5 Conclusions 

In this paper, we have studied theoretically magnetotrans¬ 
port in graphene stripes. In these stripes electrons are 
injected at one point of the boundary and focused by a 
perpendicular magnetic field onto another point of that 
boundary, see figure We have calculated by the NEGF 
method the generalized Hall resistance as a function of 
the magnetic field, see figure In weak fields equidistant 
focusing peaks appear, which correspond to classical cy¬ 
clotron orbits ( [ 2 ^ , see figure]^ When the magnetic field is 
increased, anomalous resistance oscillations are observed, 
which cannot be explained by classical cyclotron motion. 
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Figure 12. Energy resolved LDOS in the studied zigzag stripe (top row) and armchair stripe (bottom row). In the 
left column the LDOS has been averaged over the six carbon atoms of the honeycomb cells. Landau levels can be 
observed far from the edges. However, the bending of the energy bands in vicinity of the edge is seen more clearly in 
the figures of the middle and right column, where the LDOS is averaged only over the subset of atoms shown in the 
inset. In the zigzag stripe it can be seen how the zeroth Landau level (at E = 0) splits into an edge state (top, right) 
and a non-dispersive surface state (top, middle), when the edge is approached. This surface state is not present in the 
armchair stripe. 


By means of a simplified model, we have shown that 
all calculated resistance oscillations can be understood by 
the interference of the plane wave part of the occupied 
edge channels, see figure The anomalous resistance os¬ 
cillations are beatings, which appear when only some few 
edge channels are occupied and only some few plane waves 
are superimposed. Thus, the oscillations are very clear 
and distinct, if only two Landau levels are occupied. The 
frequency of the resistance oscillations increases rapidly, 
when the magnetic field is increased and a Landau level 
is depleted, because the momentum of the corresponding 
plane wave (and hence, its frequency) is also increasing 
rapidly, see figureDue to computational limitations, the 
studied graphene stripes have been relatively small and 
the magnetic field has been relatively strong. However, 
due to the good agreement of the simplified model with 
the NEGF calculations, we have used this model to show 
that our findings are expected to appear also in larger 
stripes at lower magnetic fields. As the resistance oscil¬ 
lations, classical focusing peaks as well as the beatings, 
are due to the interference of the edge channels, we have 
also given a rule of thumb (25) for the required number of 
specular reflections. 


Studying the effect of the edge shape of the graphene 
stripes on the magnetotransport, we found that a finite 
current flows on the armchair edge, whereas the current 
vanishes on the zigzag edge, see figures andBy means 
of the simplified model, the different edge currents can be 
traced back to the fact that at an armchair edge carbon 
atoms of both sublattices appear, while at a zigzag edge 
only atoms of one sublattice are present, see the inset of 
figure 1^ We have also shown in figures and that the 
number of spatially separated edge channels in the local 
current equals the number of occupied Landau levels. The 
discrete Landau levels can be seen clearly in the LDOS 
in figure However, the bending of the Landau levels 
in vicinity of the edge as well as spatially separated edge 
channels can be hardly recognized, if the LDOS is aver¬ 
aged over the six carbon atoms of the honeycomb cells. 
They can be made visible, if the LDOS is averaged only 
over a subset of the carbon atoms. 
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